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Differential positivity on compact sets 


F. Forni 


Abstract —The paper studies differentially positive systems, 
that is, systems whose linearization along an arbitrary trajec¬ 
tory is positive. Extending the results in [7], we illustrate the 
use of differential positivity on compact forward invariant sets 
for the characterization of bistable and periodic behaviors. Ge¬ 
ometric conditions for differential positivity are provided. The 
introduction of compact sets simplifies the use of differential 
positivity in applications. 

I. Introduction 

A linear system is positive if some cone /C on the system 
state space is invariant for the dynamics [3]. Positivity 
strongly restricts the behavior of a linear system. Under mild 
conditions, the ray At) £ K, given by the Perron-Frobenius 
eigenvector v £ 1C is an attractor for the system dynamics, 
[2], [3]. This fundamental property is exploited in a number 
of applications [5], [12], [15], [14] 

Differential positivity extends linear positivity to the non¬ 
linear setting. A nonlinear system is differentially positive 
if its linearization along trajectories makes a cone (field) 
invariant [7]. Differential positive systems are a large class 
of systems encompassing monotone systems [1], [9], [18]. 
Under mild conditions, the trajectories of a differentially 
positive system converge to a one dimensional attractor, 
a relevant property for the study of bistable and periodic 
behaviors. In comparison to linear positivity, this attractor 
is not a ray, but a curve, possibly closed in the presence of 
attractive limit cycles. 

In this paper we make differential positivity readily avail¬ 
able for the analysis of nonlinear systems by deriving a 
number of geometric tools. The problem of establishing the 
differential positivity of a system is encoded into a set of 
pointwise geometric conditions to test. Differential positivity 
is then used to derive novel methods for the analysis of 
simple attractors of nonlinear systems, typically capturing 
bistable and periodic behaviors. Compactness simplifies the 
use of differential positivity in applications. In particular, the 
restriction to compact and forward invariant sets makes the 
geometric conditions for differential positivity much simpler 
to verify in practice. 

Section UI1 introduces the notion of cone fields and charac¬ 
terizes their representation as a set of inequalities. The two 
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large families of polyhedral and quadratic cone fields are 
illustrated. Section [In]recalls the basic notions of differential 
positivity. Section [TV] provides a set of geometric conditions 
for testing the differential positivity of a system. The use 
of the geometric conditions is illustrated on a cooperative 
system and on a nonlinear pendulum. The role of differen¬ 
tial positivity for the analysis of the asymptotic nonlinear 
behavior is discussed in Section 0 and illustrated on the 
Kuramoto model in Section m Conclusions follow. 

II. Cone fields 

We recall some basic geometric notions on Riemannian 
manifolds which will be useful to the discussion on differ¬ 
ential positivity. Let A" be a smooth //-dimensional manifold 
endowed with a Riemannian metric (■, ■) x : T x XxT x X -A R 
where T x X denotes the tangent space at x £ X. We will use 
TX to denote the tangent bundle of X, and \Sx\ x to denote 
\J(5x, 8x) x for all Sx £ T x X. Furthermore, given any 
function ip : X —I- y between manifolds, dip(x) : T x X -A 
T^( x )y will denote the differential of ip at x £ X. Given any 
set S C T x X we will write dip(x)S := {ip{x)5x \ Sx £ 5}. 

To extend linear positivity into the differential setting, we 
will exploit the notion of conal manifold, that is, a manifold 
X endowed with a cone field 

JC £ {x) C T x X Vx £ X. (1) 

where 0 < e <C 1. At any x, IC £ (x) is just a cone in the 
vector space T x X. The role of the parameter e is clarified 
by the property IC £2 (x) \ {0} C K £l ( x ) \ {0} if £i < e 2 , that 
we assume throughout the paper. We use AC(x) := K,q{x). 

For each e, each cone K, e (x) is closed and solid , and 
satisfies (i) AC e (x) + /C e (x) C IC £ (x ), (ii) aJC £ (x) C /C e (x) 
for all a > 0, (iii) /C £ (x) IT — /C e (x) = {0}, which make 
IC £ (x) convex and pointed. To avoid pathological cases, we 
assume that for every Xi , x 2 £ X, there exists a linear 
invertible mapping F(xi ■ x 2 ) : T Xl X -A T X2 X such that, for 
any s, T(xi,x 2 )IC £ (xi) = IC £ (x 2 ). Cone fields satisfying 
this property are said regular. 

To make use of a cone field in computations, the cone 
field K, e {x ) will be represented by 

Sx £ /C e (x)\{0} Ki fx, f X ) >£, Vi£lCN (2) 

V I &X\xJ 

where each Ki : TX —> R is a smooth function. For 
simplicity, in what follows we will use K > e to denote 
the component-wise inequality Ki > e for each Vi £ T. 
For the cone AC(x), 0 is equivalent Ki(x,Sx) > 0. The 
normalization in (|2) maps every ray A Sx £ AC(x), A > 
0, into the point of the unit sphere {v & T x X \ \v\ x = 




1}. The normalization makes the representation independent 
of the length of the tangent vectors 6x. Examples of smoothly 
varying polyhedral and quadratic cone fields are provided at 
the end of the section. A simple illustration is in Figure [l] 
A cone field carries naturally the useful notion of conal 
curve 7 : R —>■ X, which is an integral curve of the cone 
field 

j(s ) € lC(y(s )) Vs S K . (3) 

We make the standing assumption that | 7 (s)| 7 ( s j = 1 for any 
s £ R. From ([2]), a conal curve satisfies K(j(s), j(s)) > 0. 

K{x) 

) 





Fig. 1. Polyhedral and quadratic cones at T x X. 

Example 1: polyhedral cone fields. Fet T be a smooth 
manifold of dimension n endowed with a Riemannian metric 
(•, •) x . For any i £ X := {1,..., rn) C N, m > n, define 

Ki(x,6x) := (Fi(x),6x) x (4) 

where Ffix) £ T x X \ {0}. We assume that: A1 for every 
pair of points x \. x 2 £ X, there exists a linear invertible 
isometry T(x i,x 2 ) : T Xl X —> T X2 X such that FfiX 2 ) = 
T{xi,x 2 )F i {xi)\ A2 {Fi,...,F m } is a smooth full-rank 
distribution; A3 the set of constraints 0-01 is feasible for 
some e > 0. Then, for all 0 < £ < e, the cone field IC e (x) 
given by <|2),(|4) is solid, pointed, convex and regular. 

Assumption A3 guarantees that, for e < e, IC e (x) 
is solid, pointed and convex by construction, since the 
distribution is full-rank. For regularity, consider any pair 
of points xi,x 2 £ X, and take (F^x 1 ), ■ ) X1 = 

£ and Sx 2 = T{x\,x 2 )8x\. Then, (Fi(x 2 ),Sx 2 ) X2 = 
(T(x 1 ,x 2 )F i (x 1 ),T(x 1 ,x 2 )6x 1 ) X2 = (Ffix i),6xi) Xl = 
e\5 Xl \ Xl . Thus, {Fi(x 2 ),js^) X2 = = e since 

T(xi,x 2 ) is an isometry. j 

Example 2: quadratic cone fields. Fet ^ be a smooth 
manifold of dimension n endowed with a Riemannian metric 
(•, -) x . Consider m> n vector fields Ffix) £ and 

define 

Ki(x, 6x) := ( F 1 (x),Sx) x 

K 2 (x,6x) := ( F 1 (x),6x)l - (F i (x)fix) x (F j (x),5x) x . 

(5) 

Assumptions A1-A3 with the additional condition A4 
(Fi(x), Fi(x)) x = 0 for any i > 1, guarantee that the cone 
field ([ 2 ]), ([5]i is solid, pointed, convex and regular for £ < £. 

Note that K 2 > 0 characterizes a double cone which is 
refined to a pointed cone by Ki > 0. The cone fields AC and 
KL e are solid and convex by construction. Regularity follows 
from the observation that {Ffix 2 ), 5x 2 ) X2 = (Ffix 1 ), fei) Xl 
for each £ X and each 6x 2 = T{x\,x 2 )8x\. j 


III. Differential positivity 
A. Differential positivity in forward invariant regions 

A linear system x = Ax, x £ is positive if there 
exists a cone K, C R" which is forward invariant for the 
system dynamics, i.e. e At IC C K. for t > 0, [3], Differential 
positivity extends linear positivity to nonlinear dynamics 

E : x = f(x) x £ X, ( 6 ) 

by requiring that a given cone field is forward invariant for 
the prolonged dynamics [4], 

52 ix - Sx)£ Tx ■ <7) 

0 represents the linearization of E along its trajectories. 
For simplicity we assume f £C 2 and forward completeness 
of E. We use if{t,x 0 ) : R x X —> X to denote the state 
reached at time t by the trajectory of E from the initial 
condition xq. Indeed, f)(-,Xo) £ E. We also use ip t (•) : X —> 
X, which maps any x £ X into ipt(x) := fi>(t,x). Finally, 
for any (x, 8x) £ TX, note that (ijj(-,x),d x fr(-,x)6x) £ ST,. 

Revisiting the definitions in [7], consider any forward 
invariant region C C X. We say that E is differentially 
positive in C with respect to the cone field KL{x) if for all 
x £ C and t > 0, 

dif t (x)JC(x) C JC(tf t (x)) . ( 8 ) 

Differential positivity captures the invariance of the cone 
field along the linearized dynamics. Furthermore, we say 
that E is (uniformly) strictly differentially positive if it is 
differentially positive and there exists T > 0 and e > 0 such 
that, for all x £ C and t >T, 

dfit(x)K.{x) C K, e {fi> t {x)) . (9) 

Strict differential positivity captures the contraction of the 
cone field K.(x) along the linearized dynamics, as shown in 
Figure [2] 


M x ) 


IC(x) 


I df t {x)K{x) 
tC{f t (x)) 


Fig. 2. Differential positivity: forward invariance of the cone field. Strict 
differential positivity: contraction of the rays of the cone field. 


B. Contraction of the Hilbert metric 


The contraction of the cone field along trajectories has 
a metric characterization based on the Hilbert metric, [3]. 
From [7, Section VI], for any given x £ C, take any 6x, 6y £ 
K{x) \ {0} and define 5y) := inf{A > 0 | X6y — 

8x£lC(x)} PI and mic( x )(8x, 6y) := sup{A > 0 | — A 6y £ 

K.[x)}. The Hilbert metric d/q^) reads 


dfc(x)($x,6y) := log 


/ M K{x) {6x,Sy) \ 

V rn K(x){dx, 8y)) 


( 10 ) 


1 (5: r, 5y) := oo when {A E M>q | A Sy — 5x E /C(ai)} = 0. 







dic{x) measures the distance between rays of the cone. 
It defines a metric in IC(x) fl {Sx £ T x X \ \6x\ x = 1}. 
Furthermore, for any Sx,5y £ K,(x), dfz{ x ){Sx,Sy) = 0 if 
and only if Sx = A Sy with A > 0, and d/c( x ){aSx, (35y) = 
djc( x ){5x , Sy) for any a > 0 and /? > 0 . 

The contraction of the cone field along trajectories is 
captured by the exponential convergence of the Hilbert 
metric, as stated by the next lemma. 

Lemma 1: [7, Theorem 2]. Let S be a strictly differen¬ 
tially positive system with respect to the cone field K,(x) in 
the forward invariant set CCA. Then, there exist k > 1 and 
A > 0 such that, for all x £ C, Sxi, Sx 2 € K.(x), and t > T, 

dic{M^)(dtpt(x)Sx 1 ,dipt(x)8x 2 ) < ke~ x(t ~ T) A (11) 

where A := sup{d/c(ao(*h> v 2 ) | vi,v 2 £ JC e (x)} < 00 . j 

Lemma |T] and the following (mild) technical assumption 
are crucial for the theorems of Section 0 
Assumption 1: (IC(x)r\{Sx £ T x X \ \Sx\ x = 1},<7k:( x )) is 
a complete metric space for all x £ C. 

The reader is referred to [3, Section 4], [10, Section 2.5], or 
[19] for examples of complete metric spaces on cones. 


IV. Geometric conditions 

To provide geometric conditions for differential positivity 
we first reformulate the property using the representation Q. 
For instance, E is differentially positive if for all t > 0 and 

all (x, Sx) £ TX, 

K(x,Sx) > 0 => K(ift t (x), di/>t(x)5x) > 0 . (12) 


In addition, E is strictly differentially positive if there exists 
T> 0 and e>0 such that, for all t>T and all (, x, Sx) £TX, 


K(x, Sx) > 0 => K 



dj> t (x)Sx A 
| dtp t (x)Sx\ Mx) ) 


>£. (13) 


© and © capture the invariance and the contraction of 
the cone field along trajectories, which leads to the following 
pointwise geometric conditions for differential positivity. 

Theorem 1: E is differentially positive in the forward 
invariant set C C X with respect to /C(x) if for any (x, Sx) £ 
TX such that K(x, Sx) > 0 and x £ C, 


Ki(x,Sx) = 0 => dKi(x,Sx) 


fix) 

df(x)8x 


> 0 . (14) 


Proof: Whenever (ipf,x),d x ip(-,x)Sx) £ ST, reaches 
the boundary of the cone field at time t, we have 
K(i(j(t,x),d x 'ip(t,x)Sx) > 0 and iQ(t/>(f, x), d x ip(t, x)Sx = 
0, for some i £ I. From ( fl4| ), jjfKi{ r f (t,x),d x 'f (t,x)5x > 
0, thus Ki(f>(t,x),d x tl)(t,x)Sx) either grows positive or 
remains at zero. ■ 

For strict differential positivity we need to take into 
account vectors on the unit sphere. For instance, for any 
Sx £ T x X, consider d := From ([7}. 

d = (df(x) — A(x, d)) d (15) 


where A(x, d) normalizes the action of the operator df{x) to 
guarantee that any trajectory d(-) of ©from |t?( 0 )| x(O ) = 1 


satisfies |'i?(f)| x (t) = 1 for all t > 0. For example, given 
the representation \Sx\ x := (5x T G{x)5x )^, where G[x) 
is the Riemannian tensor in local coordinates, X(x,d) := 
\d T (G(x)df(x) + df(x) T G(x) + Yfi=i d xiG(x)f (x)i) d. 
(f(x)i is the ith component of the vector f(x)). 

Theorem 2: E is strictly differentially positive in the for¬ 
ward invariant set C C X with respect to KL{x) if there exist 
T > 0 and e > 0 such that, for any (x, d) £ TX that satisfies 
K(x,d) > 0, x £ C and |#| x = 1, 


0 < Ki ( x , d) < e => 


dKi(x,d) 


fix) 

(df(x)-\(x,d))d 



(16) 

J 


Proof: Consider any trajectory (^(-, x), x)Sx) £ 

ST. For simplicity, define d(t) := |g $(tx)&x \^ — anc * 
x(t) := tp(t,x). Differential positivity follows from the 
argument of Theorem [1] For strict differential positivity, 
suppose that for all t £ [0, T], 0 < Ki(x(t), d(t)) < e. Then, 
( fltl i guarantees that Ki(x(T),d(T)) = Ki(x(0),d(0)) + 
fo tpd{i{x(t),d(t))dt > f 0 ifdt = £. A contradiction. ■ 
Because of the normalization, verifying ( fl 6 | may be a 
daunting task in practice. However, those conditions become 
simpler when the forward invariant region C is compact. 

Theorem 3: E is strictly differentially positive in the com¬ 
pact and forward invariant set C C X with respect to K.(x) 
if for any (x, d) such that K(x, d) > 0, x £ C and |i?| x = 1, 


Ki (x, d) = 0 


dKi(x, d) 


fix) 

df(x)d 


> 0 . 


(17) 


Proof: (i) For K(x,d) > 0, consider K t (x,d) = 0. 
Then, Ki {x,pd) = 0 for any p > 0, thus [()$K t (&,$)]$ = 
0. Since C is a compact set, by continuity, there exists k\ > 0 
(sufficiently large) and e > 0 (sufficiently small) such that 


0 < Ki (x, d)<£ => [dtfKi (x,d)]d < ki£ . (18) 

(ii) Exploiting the compactness of C again, © guar¬ 
antees that there exist k 2 > 0 and a small e > 0 such 
that Ki(x,d) = 0 => dKi(x,d)[ f(x) T df(x) T d] T > 
k 2 + Xk ± £ where A := max xeC> |^| a!=liA: (x , 1 ?)> 0 A (a:, d) and 
A(x, d) refers to ©. By continuity, for £ sufficiently small 
there exists > 0 such that 


0 < Ki (x, d) <£ => 


dKi(x,d) 


fix) 

dfix)d 


>k 2 + \k\£ — k$£ > + Afci . 

(19) 


for any x £ C and |$| x = 1 such that IC(x, d) > 0. 

(iii) We follow now the proof of Theorem [2] combining 
(i) and (ii). Consider any trajectory (ipf, x), d x tp(-, x)5x) £ 
ST. define d(t) := , and ~ ^ x )• 
Take T := || and suppose that for all t £ [0,T], 0 < 
Ki(x{t),d{t)) < £. Combining ( fl5] l, ( fl8| ) and ( fl9| ), we get 

f t Ki(x(t),d{t)) > 

> §, + Xk\£ - X(x(t),d(t)) [d^Ki(x(t),d(t))]d(t) > f 

s ---^---^ 

< a <feie (20) 















which leads to a contradiction by integration over the interval 
[0,T], as in the proof of Theorem [2] ■ 

Theorem [3] is illustrated in the following examples. 

Example 3: Cooperative systems. A nonlinear system £ 
given by x = f(x), x £ R", is cooperative if its Ja¬ 
cobian df(x) has nonnegative off-diagonal elements, [18]. 
Cooperative systems are differentially positive systems in 
X := R" with respect to the constant polyhedral cone 
field given by the positive orthant, that is, KL{x) := R" for 
any x £ X (Theorem [lj. Given any compact and forward 
invariant region C C R", £ is strictly differentially positive 
if the off-diagonal elements of df(x) are strictly positive 
(Theorem[3]l. To see this, consider the standard inner product 
in R” and the cone field 0 given by Fi(x) := e, for each 
x £ C and i £ X := {1,... ,n}, where e* is the canonical 
base. Then, 5x £ K,(x ) reads ef 5x > 0 for all i. o reads 
ef Sx = 0 => efdf{x)Sx > 0 for \ 8 x \2 = 1, which is 
equivalent to efdf(x)ej > 0 for all i. j £ I, j. _i 

Example 4 : Differential positivity of the pendulum. Con¬ 
sider the nonlinear pendulum given by the equations d = v, 
i) = — sin(t?) — kv + u. The linearization reads Sd = Sv, 
Sv = — cos {d)5d — k 8 v. Theorem [3] guarantees that for any 
k > 2, any input u, and any compact and forward invariant 
region C, the pendulum is strictly differentially positive in C 
with respect to the cone field 0 given by 

Sd > 0 Sd + Sv > 0 (21) 

For instance, for | [w 8 v] T | =1, (i) 8 d = 0, Sd + Sv > 0. 
Then Sd = Sv > —Sd > 0. (ii) Sd > 0, Sd + Sv = 0. Then, 
Sd + Sv = Sv — cos(d)5d — kSv > (k — 1 — l)<5t? > 0. 

Revisiting [7, Section VIII], Theorem[6]below can be used 
to establish the existence of limit cycles for u > 1. j 

V. Asymptotic behavior 
A. Stable attractors 

It is well known that under mild conditions almost every 
bounded trajectory of a monotone system converges to a 
fixed point [8], [18]. This fundamental result has been 
recently revisited in [7, Corollary 5], with a new proof based 
on differential positivity. The next theorem extends this result 
to a larger class of systems. 

Given a compact set C C X we say that a conal curve 
7 : R —> X intersects the boundary of C twice if for any 
s £ R such that 7(s) £ C, there exists an interval s<s<s 
such that 7(s),7(s) C. 

Theorem 4 : Under Assumption [T] consider a strictly dif¬ 
ferentially positive system £ with respect to the cone field 
K[x) in a compact forward invariant region C C X. Suppose 
that every conal curve 7 : R —> X intersects the boundary of 
C twice. Then, from almost every initial condition in C the 
trajectories of £ converge asymptotically to a fixed point, j 

The next theorem exploits the combination of differential 
positivity with the existence of an invariant vector field 
v{x) £ lC e (x). It shows that the trajectories of the system 
converge asymptotically to a one dimensional attractor given 
by the image of an integral curve of v{x). 


Theorem 5: Under Assumption [I] consider a strictly dif¬ 
ferentially positive system £ with respect to the cone field 
K[x) in a compact forward invariant region C C X. Let 
e > 0 and suppose that there exists a complete vector field 

v(x) £ IC e (x) \ {0} such that 


limsup \df> t {x)v{x)\^ x ) < 00 ; (22) 

t—f OO 


v(tpt{x)) 


dif t {x)v{ x) 
\di> t {x)v(x)\^ x ) 


\/x £ C,\/t > 0 . (23) 


Then, there exists an integral curve of v(x) whose image is 
an attractor for all the trajectories of £ from C. j 

Finding the vector field that satisfies ( |22| ) and ( [23] ) can 
be difficult in general. However, the presence of symmetries 
in the system makes ( [22] ) and ( |23| ) tractable conditions. An 
example is given by consensus dynamics [13], [16] where 
n agents communicate with their neighborhoods to achieve 
consensus, typically given by the manifold x 1 = • • • = 
x n . The invariance of the consensus manifold dictates the 
symmetry df(x) 1 = 0, 1 := [1... i] T £ T x X, which 
makes ([22] and ( |23| trivially verified by v(x) := J-. This 
observation is used in Section [Vi] for the analysis of the 
Kuramoto model. 

Replacing v(x) in Theorem [5] with the system vector field 
f(x), the next theorem gives conditions for the existence 
of attractive limit cycles. This result is compatible with 
Theorem [5] since conal curves can be closed curves. 

Theorem 6 : Under Assumption [I] consider a strictly dif¬ 
ferentially positive system £ with respect to the cone field 
IC(x) in a compact forward invariant region C C X. Suppose 
that C does not contain any fixed point. If f{x) £ A3 e (tc)\{0} 
for any x £ C, for some s > 0, then there exists a unique 
attractive periodic orbit contained in C. j 

Remark 1: Converse results for hyperbolic limit cycles 
can be found in [11]. Theorem [6] revisits [7, Corollary 2], 
which requires differential positivity in the whole manifold 
X, a condition weakened by Theorem [6] The key step for 
this result is a new proof that does not use the so-called 
Perron-Frobenius vector field of [7, Section VI]. _i 


B. Proofs 

Proof of Theorem [7] Suppose that for some x £ X, the 
trajectory if{-,x) does not converge to a fixed point and 
denote by w( x) the w-limit set of x. Then, there exists a 
sequence of time instant tk —> 00 as k —» 00 and c > 0 such 
that \f(if(tk,x))\ > c, for all k £ N. By continuity, since 
ip(t, x) £ C for all t > 0, there exists a small constant p > 0 
such that \f(xf{t,x))\ > c for all t £ [tk — p,tk + p\- 

Recall also that (f)f,x), f(rf)(-,x)) £ <5£ since 

jtf{ip(t,x)) = df{if(t,x))f(tp(t,x)). Thus, either (i) there 
exists r > 0 such that f(if(t,x)) £ —K.{ip{t,x)) U 
IC(fj(t,x)) for all t > t, or (ii) f(f>(t,x)) ^ -~K,(if(t,x)) U 
x)) for all t > 0. 

For (i), note that the points ipt+ T {%) for t > 0 belong to the 
image of a conal curve, since ^fip t+T (x) = /('0 t+T (x)) £ 
Kfift+rix)) (or -f(ip t+T (x)) £ ICfift+rix))) for any t > 0. 
Then, since \f(i/j t + T (x))\ > c for all t + r £ [tfc — p,tk + 



p\, exploiting the fact that every conal curve intersects the 
boundary of C twice, there exists a time T > 0 such that 
f. C. This contradicts the forward invariance of C. 

For (ii), we show that the basin of attraction of 
A := oj( x) has dimension n — 1 at most. We need 
a preliminary result. Consider any Sx £ /C e (x) for 
some e > 0. Then, for a > 0 sufficiently large 

/( x) + aS:x £ /C(x). Thus, by projective contrac¬ 
tion, limd K ^ t(x ))(dip t (x)[f(x) + aS x\,dif t (x)aSx) = 
lim d K f^ t r x)) (adif t {x)8x + f(if t {x)),adif t {x)5x) = 0. 
Since \ffft k (. x ))\ > c, it follows that 

^lim \dip t (x)5x\ 1 p t ^ x ' ) = oo \/Sx £ int IC(x) . (24) 

By contradiction, suppose now that the basin of attraction 
B a has dimension n. Consider a conal curve 7 such that 
7(0) £ A and 7 (s) £ S 4 \ {A} for all s £ [ 0 ,s] C 
R. By assumption, for all s £ [0,s], ift(y(s)) converges 
asymptotically to A as t — > 00. Moreover, f{ift(y{s))) ^ 
1 Cff t (p(s))) for all t > 0, and \fff tk (t(s)))| > c for k 
sufficiently large. From ( |24| ), it follows that 

lim \dif t {x)y(s)\ Mlis)) = 00 • ( 25 ) 

£—>00 v v 

From ( [25] ), the length of the curve if t { 7(-)) grows unbounded 
as t ->• 00, since Ifftfjf))) = In I £M7(s))ds\ MyM ) = 
f 0 \d'i/jt.('y(s))j(s)\^ t ( 7 ( s ))ds. It follows that there exists 
T > 0 such that ifrilis)) ^ C, contradicting the forward 
invariance of the set. ■ 

Proof of Theorem [j] [(i) boundedness] Under the assump¬ 
tions of the theorem, we prove that for any given Sx £ T x X, 


lim \dif t (x)Sx\^ x ) < 00 . (26) 

Consider first the case dift(x)5x —JC(ift(x)) U 

K.(ift(x)) for all t > 0. Then for a > 0 suf¬ 
ficiently large 8x + av(x) £ KL(x), thus Lemma [T] 

gives \\m d K t^ t t x )){df) t (x)[av{x) + 8x\,df) t {x)v{x)) = 
£—>00 

lim d K ,ij, t i x \ ) (adif t (x)v{x)+dift(x)5x,adift{x)v{x)) = 0. 
£—>00 

Since \di/jt{x)v(x)\^ t ^ x ' ) is bounded, it follows that 

lim \dif t (x)8x\ifi t r x \ = 0 . (27) 

£—>00 

Consider now the case dipT{x)Sx £ K,(iPt(x)) for some 
T > 0. If dipT{x)8x £ —K.(iI>t(x)), consider the sym¬ 
metric trajectory dipT(x)[—Sx] £ JC(i/}t(x)). Suppose that 
^lim \d'f>t{x)5x\^, t ^ x ) = 00 . Then, by projective contraction 

dK^ t+T (x)){dfit+T(x)Sx,dft+T{x)v{x)) = 0. Therefore, 
for t sufficiently large, dipt+T{x)v(x) is almost parallel 
to df> t +T{x)8x. By linearity of dip t +T(x), it follows that 
dipt+T(x)v(x) grows unbounded, contradicting ( [22] ). 

[(ii) horizontal contraction] We show that any pair of 
points in C converge to the image of an integral curve of 
v(x). Consider any curve 7 o(-) : [0,1] — > C connecting two 
different points 70 ( 0 ), 7 o(l) £ C. The evolution of the curve 
along the flow is given by 7 t (s) := Note that 

£lt(s) = aV> t (7(s))7(s). From (i), either lim ^7 t (s) = 0 

by eg, or there exists T such that ^7 t (s) £ — /C(7t(s)) U 


/C( 7 t (s)) for all t > T. In this last case, using ^26) , 
is bounded, and d K{nt(s)) {-^'Yt{s),dfj t ( r y(s))v('y(s))) = 0. 
Thus, in the limit of t —► 00 , ^ 7 t (s) becomes parallel to 
u( 7 t(s)). It follows that f’ti'yis)) belongs to the image of an 
integral curve of v(x). 

[(iii) uniqueness] By contradiction, let A\ and A 2 be 
images of two distinct attractive integral curves of v(x). Take 
X\ £ A \, x 2 £ A 2 , and consider a new curve 7 connecting 
them. Along the flow, iptilis)) converges to an integral curve 
of v(x). Thus, by completeness, Ai and A 2 are subsets of 
the image of the same integral curve of v(x), contradicting 
the initial hypotesis. ■ 

Proof of Theorem [ 6 ] From the conditions of the theorem, 
there exists 0 < C\ < C 2 such that c\ < |/(x)| < C 2 for 
all x £ C. Reasoning like in the proof of Theorem 0 by 
replacing v(x) by /( x), we get 

lim |dV , t(a ; )<k 4 l />t(x) < 00 (28) 

£—>■00 

for any x £ C and Sx £ T x X. Furthermore, we can conclude 
that there exists a unique integral curve of the vector field 
f(x) - a trajectory - whose image is an attractor for every 
trajectory of E from C. In what follows we show that this 
curve is closed, by following the proof of [7, Corollary 2], 

Let iff, x) be such a curve, for some x £ C. iff, x) does 
not converge to a fixed point and belongs to a compact set, 
therefore there is a point x* = if ft* ,x) whose neighborhood 
B p (x*) is visited by the trajectory infinitely many times, for 
any p > 0. Take p > 0 sufficiently small and let S C C 
be a transversal manifold to the trajectories of E in B p (x*), 
namely ffx) £ T X S for all x £ S D B p (x*) and x* £ S. 
Consider a sequence of time instants t-k —t 00 for k -£ oc 
such that if(tk,x*) £ S (1 B p (x*), and the subsequence of 
time instants tj. —> oc for j —> oc such that if(fk ,x*) £ 
SHBe{x*), 

Consider now any curve 7 (-) : [0,1] —> S n B p (x*) 
such that the length of ifyf)) = f 0 | 7 (s)| 7 ( s )rfs is less 
than or equal to p £ I and 7 ( 0 ) = x*. By ( |28| ), there 
exists C 3 > 1 such that limsup^^ £(ift(jf))) < C 3 P- 
Furthermore, lim^^ d K(i , t(x ^(£if t ('y(s)), ffft{l{s))) = 
0 for all s £ [0, 1 ], that is, in the limit of t —t 00, the 
image of iftf)f)) converges asymptotically to the image 
of the attractive integral curve of the vector field ffx). 
Thus, because of the transversality of <S, the combination 
of the bound limsup^^ £fftflf))) < c$P and of the 
convergence of the Hilbert metric guarantees that there exist 
k sufficiently large and 77 (s) typically small, such that 
if{tk + T k (s),y(s)) £ S n Be(if(t k ,x*)), Vs € [0,1]. 
Thus, there exists a j sufficiently large such that ifftk + 
rfe (s), 7 (s)) e <Sn2?p (t/’fifc, C SnB 2 p (x*). It follows 
that the return map on S is necessarily a contraction, which 
implies that the attractor is a closed curve. ■ 


VI. Extended example: all-to-all Kuramoto 
We study the synchronization of Kuramoto dynamics for 
the case of all-to-all coupling with homogeneous velocities 
at zero using the results of Section [V] Differential positivity 


is verified by using the geometric conditions of Section IV 



Consider the interconnection of n agents (phase variables) 
given by i9 k = i Yf=i sin (^ - $fc) where d k £ S. The 
synchronization manifold is given by A := {•& £ S n | = 

• • • = $„}, where $ := [tfi ... ] T , k £ X := {1,..., n}. 

Using 1 := [1 ... i] T , the prolonged dynamics reads 


balanced phase (p = 0) or saddle point (Efcgz s i n (^fc — 
(ft ) 2 = 0 ), there exists a sufficiently small e > 0 such that 
sin(r?fc — </>) 2 > e and p > 0 for every i9 £ C. Thus, 
there exists a sufficiently large A such that <W) = 0 => 

&K 2 (0,60)> 0 . 


r 0 = ±S(0)1 
\ M = £C(i?)M 


(0, <h?) g TS” (29) 


where for k,i £ X, S k i(fi) := sin($i — $&), Cfcfc($) := 
- COS and C ki (i9) := cos(i?j-i?fe) for fc ^ j. 

The invariance of the synchronization manifold is captured 
by the identity C(0)1 = 0. The constant vector field i>(0) = 
- 7 = satisfies the conditions of Theorem |5j Thus, exploiting 
the theorem, given a compact and forward invariant region 
C, all trajectories from C converge asymptotically to the 
synchronization manifold A C C, provided that the system 
is strictly differentially positive with respect to some cone 
field /C($) such that 1 £ int /C(i9) (interior of /C(0)). 


Define II := 


I 



and consider the centroid p (. ;i0 := 


n Efcex e ™ k ■ P £ [0,1] is a measure of the synchrony of the 
phase variables [17, Section III]. We show that the system 
is strictly differentially positive in any compact and forward 
invariant region C C S n with respect to the quadratic cone 
field ([5) given by 


JTiOMtf) := 1 T <W 

K 2 (i9, Si?) := e 2 Ap Si9 T ll T Si9 - Si9 T nSi9 


where A > 0 is a parameter defined below. An illustration is 
provided in Figure [3] 



Fig. 3. Increasing p widens the cone. The cone is widest at p = 1. 

For K 2 ( 1 ?, S'd) = 0, using the identities l 1 C(r?) = 0 
and IIC($) = C(i?), S’d) along the flow of the 

prolonged dynamics reads 

2\pe 2 Xp 8d T H T 5d - <W T (C(r?) T + C(t9))St9 (31) 

where, following [17, Section III] and [6, Section VII], 
p = p Efcgi s i n (^fe — 4 >) 2 - In particular, p = 0 for 
p = 0 (balanced phases, max spread on the circle) or for 
Efc e iSin(^fe — <f> ) 2 = 0, which occurs on isolated critical 
points given by n — to phases synchronized at <f> + 2jn and 
to. phases synchronized at <j> + tt + 2 jn, for j £ N and 
0 < to < ^ (for an extended analysis see [17, Section III]). 

For £ S” /2 := {t? £ S n | \i) k - -dj\ < f , Vfc, j £ X }, 
the quantity —S’d T (C('d) T + C(r?))<5$ > 0. For •& £ 
we can design A so that the first term 2Ape 2Ap <5i? T ll^ 5$ in 
© dominates —^ T (C(^) T H-CGiven any compact 
and forward invariant set C C § n that does not contain any 


VII. Conclusion 

Differential positivity provides a number of methods for 
the analysis of simple attractors of nonlinear (closed) sys¬ 
tems, capturing bistable and periodic behaviors. The intro¬ 
duction of compact sets simplifies the use of differential 
positivity in applications. The theory have been illustrated 
on several examples based on polyhedral and quadratic 
cone fields. Future research directions will study differential 
positivity on open systems. 
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